Neutrino Lump Fluid in Growing Neutrino Quintessence 
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Growing neutrino quintessence addresses the why now problem ol dark energy by assuming that 
the neutrinos are coupled to the dark energy scalar field. The coupling mediates an attractive force 
between the neutrinos leading to the formation of large neutrino lumps. This work proposes an 
effective, simplified description of the subsequent cosmological dynamics. We treat neutrino lumps 
as effective particles and investigate their properties and mutual interactions. The neutrino lump 
fluid behaves as cold dark matter coupled to dark energy. The methods developed here may find 
wider applications for fluids of composite objects. 



I. INTRODUCTION 

The observed accelerated expansion of the Universe 
can be described by a dark energy component [3, Q • Its 
energy density dominates that of matter at present, while 
it constituted a very small fraction of the energy budget 
in earlier stages of the cosmic evolution This "why 

now" problem has motivated the idea of dark energy be- 
ing dynamically coupled to other cosmological species. 
It has been proposed that a dependence of the neutrino 
mass on the dark energy scalar field, the cosmon, may 
naturally trigger the onset of accelerated expansion in re- 
cent times [a, la] . The background evolution of the result- 
ing cosmological model, growing neutrino quintessence, 
is similar to the concordance model with a cosmological 
constant A. 

Since the energy density in neutrinos is small, the 
cosmon-mediated attraction between neutrinos has to be 
substantially stronger than the gravitational one in order 
to be effective. This results in a fast formation of neutrino 
lumps of the size of clusters or larger at redshift around 
one. The dynamics of the perturbations in the cou- 
pled cosmon-neutrino fluid is complicated. In contrast 
to models of uncoupled or weakly coupled dark energy, a 
mere analysis of the background equations together with 
linear perturbation theory is insufficient. Linear pertur- 
bation theory breaks down even at large scales Q, and 
the nonlinear evolution exerts significant backreaction ef- 
fects on the background evolution. This has lead to the 
development of a specifically designed N-body based sim- 
ulation method, which accounts for local cosmon pertur- 
bations, relativistic neutrino motion, and backreaction 
effects [8| . These simulations are, so far, successful until 
z f=3 1, where a collection of spherical neutrino structures 
has formed, cf. Fig. [TJ 

Although it is numerically challenging to resolve the 
internal dynamics of the neutrino lumps, these details 
may not be crucial for the broad cosmological picture. 
In gravity, e.g., the detailed evolution inside galaxies or 
clusters is not relevant for the cosmological evolution. 
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FIG. 1. Neutrino structures in a simulation box of comoving 
size L — 600/i _1 Mpc at redshift 0=1. Shown are regions 
with a neutrino number density contrast above 5 @]. 



Once the neutrino lumps have formed, one would like to 
use a picture of a pressureless fluid of neutrino lumps. 

In contrast to the universal properties of gravity, where 
only the total mass of a bound object matters, the under- 
standing of neutrino lumps needs more information. The 
mass of a lump with a given number of neutrinos is still 
expected to depend on the local value of the cosmon field 
(p averaged in a region around the lump. This effective 
coupling of tp to the lumps induces an effective attrac- 
tive interaction between the lumps. Since the lumps are 
highly nonlinear objects, the (^-dependence of the mass 
is sensitive to the total number of neutrinos in the lump 
and possibly even to additional properties of the lump. 

This work presents analytical and numerical studies of 
the properties of the neutrino lumps. We indeed find an 
effective description. This opens the possibility for an 
approximate and much simpler approach to the under- 
standing of the cosmological evolution of growing neu- 
trino quintessence for the period after the formation of 
the lumps. 

The paper is organized as follows. We collect some 
basics of growing neutrino quintessence in Sec. |TT] and 
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motivate the approach taken in this work. Section IIIII 
describes the effective cosmological dynamics in the pres- 
ence of stable neutrino lumps. Starting from the basic 
idea of approximating lumps as particles, we eventually 
develop a simplified simulation scheme of growing neu- 
trino quintessence. Some more technical aspects required 
for this scheme are postponed to Sec. IIVI The question 
of stability of neutrino lumps is discussed in Sec. [V] We 
conclude in Sec. [VlJ 

II. FUNDAMENTALS AND MOTIVATION 

After briefly summarizing the basics of growing neu- 
trino quintessence in Sec. Ill A[ we explain the main idea 
of this work. We give physical arguments why the neu- 
trino lumps may be approximated as nonrelativistic par- 
ticles. This forms the basis of the effective description of 
the cosmological dynamics presented in Sec. IIIII 

A. Basics of growing neutrino quintessence 

The cosmon-neutrino coupling is described by the 
energy-momentum exchange 

V A T$ = +/?T W 9'V, (1) 

v A rg = -pr^Wif, (2) 

where /3 is a dimensionless coupling parameter and Tt v \ = 
T^g^x is the trace of the neutrino energy-momentum 
tensor. We work in units where 8ttG = 1 and use the met- 
ric convention ds 2 = -(1 + 2^)dt 2 + a 2 (l - 2§)dx 2 . This 
type of coupling corresponds to early proposals of cou- 
pled quintessence 0, [l(| ■ On the particle physics level, 
the coupling is realized as a dependence of the (average) 
neutrino mass m„ on the cosmon field 6] : 



For simplicity, we consider the case of a constant coupling 
parameter f3 as used in, e.g., @, H, [Hj]. Typical values 
are of order j3 ~ — 10 2 . 

When the cosmon rolls down its potential towards 
larger values, a negative /3 implies a growing neutrino 
mass. As long as the neutrinos are highly relativistic 
(w u ~ 1/3), the trace Tt v \ — —p v (l — 3w u ) is close to 
zero and hence the coupling is small. This changes once 
the neutrinos become nonrelativistic. The coupling then 
stops the further evolution of the cosmon resulting in an 
effective cosmological constant. In this way, the model 
addresses the "why now" problem of dark energy. As in 
standard quintessence models [12, [l]| , the energy density 
of the dark energy scalar field ip decays similarly to the 
other species during most of the cosmological evolution, 
thereby alleviating the fine-tuning of the present amount 
of dark energy. 



The energy-momentum exchange, Eqs. (fT]) and ([2]), im- 
plies [3], in the Newtonian limit, an attractive force be- 
tween the neutrinos of order 

|F|«|/?V^|«2/3 2 |F gravity |. (4) 

We shall see that the interaction between neutrino lumps 
is similar but with an effective coupling weaker than (3. 



B. Lumps as nonrelativistic particles 

The simulations of growing neutrino quintessence have 
shown that, after a phase of rapid neutrino clustering, al- 
most all cosmic neutrinos are bound in roughly spherical 
lumps, cf. Fig. [TJ 

Inside these lumps, the neutrinos have relativistic ve- 
locities [1, [nj. For the neutrino fluid alone, one thus 
observes a large pressure such that a nonrelativistic treat- 
ment is not applicable. This is reflected in the equation 
of state w v — Pv/ Pv, which reaches uv ~ 0.1 at z = 1 [8J. 
Nevertheless, we argue that the lumps as static bound 
objects behave as particles with vanishing internal pres- 
sure. The pressure induced by the neutrino motions is 
cancelled by a corresponding negative pressure of the lo- 
cal cosmon perturbations. Furthermore, the peculiar ve- 
locities of the lumps are nonrelativistic. This is similar to 
a gas of atoms at low velocities. Although the electrons 
move at high velocities, their contribution to the pressure 
is cancelled by a contribution from the electromagnetic 
field. 

Whereas the total pressure of a lump vanishes, the 
contributions of neutrinos and the cosmon perturbation 
do not cancel locally. The neutrinos are rather concen- 
trated and hence their pressure contribution is restricted 
to a small radius. The cosmon perturbation, in contrast, 
extends to larger distances, analogously to the gravita- 
tional potential around a massive object. The cancella- 
tion thus only refers to the integrated contributions at a 
sufficiently large distance from the lump. 

In the following, we discuss this in more detail. Since 
gravity is subdominant compared to the fifth force, cf. 
Eq. ((4]), it may be neglected for a simple discussion. On 
general grounds, one can show that a bound object has 
vanishing pressure if three conditions are met: 

1. The object is described by a conserved energy- 
momentum tensor. 

2. The energy- momentum tensor vanishes outside a 
volume surrounding the object. 

3. The energy-momentum tensor is static. 

The argument is given in Sec. IIV Al For the purpose of 
illustration, we have numerically simulated an exemplary 
spherical neutrino lump satisfying these idealized condi- 
tions. The staticity of the lump was realized by a hy- 
drodynamic balance equation, cf. Sec. lVBl The neutrino 
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pressure integrated to a comoving radius r from the cen- 
ter is given by a sum over particles p 

Pu{r) = J\^dr^)\Tl v)i = £ ^m p v 2 p , (5) 

with the Lorentz factor 7 p and the determinant of the 
spatial metric \/p( 3 ) w a 3 . The contribution of the cos- 
mon perturbation is 



6a 2 



(G) 



where we have subtracted the pressure induced by the 
background field (p. Figure [2] shows the cancellation of 
the pressure contributions for large radii. As already ex- 
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FIG. 2. Integrated pressure contributions P v (black dashed), 
Ps v (black dotted), and their sum (red solid), normalized by 
the lump mass Afj. 

plained, the cosmon contribution is more extended than 
the neutrino contribution. 

In the cosmological context, the aforementioned con- 
ditions are met, at best, approximately and realistic neu- 
trino lumps will not be exactly pressureless. We shall now 
discuss the three conditions. First, only the total energy- 
momentum tensor of neutrinos, local cosmon perturba- 
tion, and background cosmon is conserved. The back- 
ground field, however, cannot be attributed to the lump 
(otherwise, the second condition would not be satisfied). 
The energy-momentum tensor of the lump, defined to in- 
clude the neutrinos and the local cosmon perturbation, 
is thus not exactly conserved due to exchange between 
the lump and the outside cosmon field. Finally, even for 
a virialized lump with a fixed number of neutrinos, the 
energy-momentum tensor is not static. Due to the time 
evolution of the outside cosmon field, the mass of the 
neutrinos and therefore the mass of the lump changes. 

One may argue that these effects are suppressed by the 
difference in the relevant time scales for the dynamics of 
the lump and the cosmological evolution. Indeed, the vi- 
olations of staticity and energy-momentum conservation 



are proportional to the time derivative of the cosmon field 
averaged on length scales much larger than the size of the 
lump. This is suppressed by the fact that the associated 
time scale is large as compared to the dynamical time 
scale of the lump. The effective description of growing 
neutrino quintessence presented in the next section as- 
sumes that the pressure of neutrino lumps approximately 
vanishes. 



III. EFFECTIVE DYNAMICS 

The approach of this section is to treat the neutrino 
lumps as effective particles. We then merely have to char- 
acterize their mutual interactions and their influence on 
the background as well as on the gravitational potential. 
A numerical treatment of the internal structure of the 
lumps will no longer be required. 

In Sec. IIII A[ we shall describe how lumps can be 
treated as particles with an effective coupling. Sec- 
tion IIII Bl derives the equation of motion for these ef- 
fective particles and explains how to calculate the rele- 
vant potentials: the large-scale cosmon (p and the grav- 
itational potential Finally, we explain in Sec. IIII CI 
how the results can be used to construct the simplified 
simulation scheme for growing neutrino quintessence. 



A. Description of lumps 

Let us introduce a comoving length scale A, which is 
larger than the typical lump sizes but smaller than their 
typical distances (the mean distance between neighboring 
lumps is of order lOO/i Mpc). On scales larger than 
A, a lump I at comoving coordinates xi looks effectively 
point-shaped, 




x - Xi), 



(7) 



the amplitude A^ v being given by the integrated local 
energy-momentum tensor of the lump, 



A>™ = / ^v^CW- 



(8) 



We will see in Sec. IIV Al that this indeed reduces to the 
standard one-particle case 



(9) 



where M; is the lump's rest mass (consisting of a neutrino 
and a cosmon contribution) and is its four-velocity. 
The Lorentz factor is defined as 7 = yf—goo u . In the 
background metric, we have 7 = u°. The result for A^ v 
is a consequence of the approximate pressure cancellation 
discussed in Sec. Ill Bl 

The interactions between the lumps are mediated by 
the cosmon field tp. Given that the distances between 
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the lumps are greater than A, it suffices to consider the 
smoothed field (indicated by a hat) 



<p(x) 



■y)v(w) (io) 



with a suitable window W\ of size A. 

Analogously to the fundamental coupling parameter /3, 
Eq. (J3J, we may define the effective coupling by 



din Mi 

dip 



(11) 



The effective coupling may depend on the scale A over 
which the field is averaged. Whereas the fundamental 
coupling /3 describes the dependence of the microscopic 
neutrino mass m„ on the local cosmon field tp, the ef- 
fective coupling (3i measures the mass dependence of the 
total lump mass Mi on the large-scale cosmon value (p. 
As the fundamental parameter /3 quantifies the force be- 
tween neutrinos, cf. Eq. the effective parameter j3i 
will determine the interactions between lumps. 

We next show quantitative results for the distribution 
of lumps and the effective couplings at z = 1. For this 
purpose, we have performed 10 simulation runs with the 
method and the parameters of Ref . 0] : fundamental cou- 
pling j3 — —52, box size L = 600/i _1 Mpc, but with 
reduced resolution iV C eiis = 128 3 . The positions of the 
lumps have been identified as local maxima of the neu- 
trino density field (cf. DENMAX halo finding QJ]). A 
glance at Fig. [T] shows that there is not much ambiguity 
in identifying lumps. 

Once a stable lump has formed, the number of bound 
neutrinos is approximately fixed (neglecting merging pro- 
cesses) . It is thus natural to characterize different lumps 
by their amount of neutrinos. 

We measure the effective couplings fli and the lump 
masses Mi. The latter include a (dominant) neutrino 

contribution Mj and a somewhat smaller cosmon part 



M/ v) . 
yields 



Integration over the comoving lump volume Vi 



li M, 



with p v = \ 




(12) 



(13) 
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FIG. 3. Lump abundances, effective couplings, and lump 
masses as functions of the neutrino number fraction /; (num- 
ber of neutrinos in the lump normalized to the number of 



+ V(ip). The Lorentz factor 7 neutrinos in the Hubble volume Vh = H 3 ) at redshift 



depends on the velocities of the lumps or particles, re- 
spectively. The smoothed field p is considered external 
to the lump and thus subtracted. 

Figure [3] shows the abundance of lumps and the distri- 
butions of Pi and Mi. The couplings f3i are measured by 
numerical differentiation according to Eq. The two 

lower figures show approximate functional dependences 
on the neutrino amount with only relatively small statis- 
tical fluctuations. The effective coupling is systematically 
weaker than the fundamental coupling. This becomes 
more pronounced with increasing neutrino number. 



The error bars indicate the variance of lumps in the same bin. 



For the averaging scale A, we have taken A = 
30 /i _1 Mpc. This is clearly smaller than the typical lump 
distances ~ lOO/i^Mpc but larger than the neutrino 
concentration of the lumps. Concerning the cosmon field, 
there remains some ambiguity since we attribute only 
the cosmon perturbations at scales smaller than A to the 
lumps. If A is chosen larger, the pressure cancellation and 
thus the particle approximation are better, cf. Fig. [21 but 
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there may arise overlaps between spatially close lumps. 



B. Evolution equations 

Within our effective description, the equation of mo- 
tion of a neutrino lump is derived from the standard one- 
particle action 



-9 T i 



drMi(ip) 



(14) 



with the proper time r and the smoothed cosmon field (p 
evaluated at the lump trajectory. Along the same lines 
as for the single neutrino case [8|, we arrive at 



dvt 1 ,, „ „ 



018*$ + 0; u x d^u". 



(15) 



The left-hand side describes gravity (expansion and grav- 
itational potential), the right-hand side is due to the 
cosmon- neutrino interaction. The (spatial) term 0i^<p 
is the cosmon-mediated fifth force analogous to Newto- 
nian gravity, cf. Eq. Q. The second contribution on the 
right-hand side reflects momentum conservation: A lump 
is accelerated when it moves towards a direction where 
it loses mass. 

In order to use this effective equation of motion, we 
need to know the smoothed cosmon field ip, the gravita- 
tional potential and the background evolution. In the 
following, we shall describe how this is achieved. 

For the calculation of <p, we recall the coupled Klein- 
Gordon equation, separated in background and pertur- 
bation parts [H, 



$ + 3H<p + V'{tp) 
AStp - a 2 V"{(p)8ip 



—0T( V ), 
0a 2 5T {v) . 



(16) 
(17) 



In Eq. (TIT)) we have neglected the gravitational potential 
against the cosmon perturbation. The second equation 
is similar to the gravitational Poisson equation. A natu- 
ral choice for the cosmon potential V is the exponential 
potential V(tp) cx exp(— cap) fl6| . 

Next, we will smooth the perturbation equation (|17p . 
For the left-hand side, it is straightforward to show by 
partial integration that 

A6^p(x) = / d 3 y W\(x - y)A y 5<p(y) 



A, 



d 3 y vV^ W x {x - y) 5ip{y) 



A6<p(x), 



(18) 



up to surface terms and neglecting the metric pertur- 



bations, 



ex- 



tensor of neutrinos 



On the right-hand side, we write 
■Tt v \ with the smoothed energy-momentum 



f {v) {x) = J d 3 v^)W x {x-y)T (v) {y). 



(19) 



We next employ the relation (shown in Sec. IIV Bl) 

lumps I 

Here, the smoothed trace of the energy-momentum ten- 
sor of a lump Ti can be calculated from the effective lump 
energy-momentum tensor, Eqs. ((7|) and (|9]), T; = Tj 11 g^: 



ft = I ^yV9^Wx(x-y)Ti(y) 



li 



-W\{x - xi). 



(21) 



With these results, the smoothed perturbation equation 
eventually reads 



AS<p(x) - a 2 V"{(p)5tp(x) 



a 2 Y Pi~W\(x — xi) — 0a 2 T( v y 



(22) 



lumps I 

Assuming that all neutrinos are bound in lumps, one has 

1 



" lumps I 



(23) 



in some cosmological volume Vph ys ■ 

On scales larger than A, the window W\(x — x{) in 
Eq. (|2"2"j) may be replaced by a point cx S^(x — xi). For 
an approximate solution of Eq. (1221) at distances larger 
than A from the sources, we thus use a sum of Yukawa 
potentials, 



Sip 



_ ^ , V 47ra a; - x\\ 

lumps ( 



■am^ \x—xi 



5(f r . 



(24) 



with the scalar mass m 2 = V"(ip). The residual term 
Sip T es,i{ x ~ x i) is needed to cancel the background part 
cx T(v) on the right-hand side and to ensure Sip = in a 
simulation volume, similar to ^ros,; below. 

If the lumps are moving rather slowly compared to the 
speed of light, 72 ss 1, the two smoothed metric potentials 
are equivalent, $ f=a \F Numerically, this relation is veri- 
fied on large scales [8j . Then, we write for the smoothed 
gravitational potential induced by lumps (with the same 
approximations as for ip) 



M, 



5^(x-xi) 



Mi 



lumps I \ V 9 

The solution is (up to a constant) 



V, 



phys 



E 

lumps I 



Mi 



V 8na \x — xA 



(25) 



(26) 



where the residual contribution can be given explicitly as 
'I' res,; = Mi a2 \ x — x i\ 2 /Vphys- The total gravitational po- 
tential also includes the matter-induced potential which 
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is calculated as usual. Taking into account relativistic 
corrections would require the calculation of both poten- 
tials, * and $, cf. Ref. 0. 

In order to have a full description of the cosmological 
dynamics, we still need to describe the evolution of the 
cosmological background, i. e. the Hubble expansion H 
and the background cosmon (p. The background evolu- 
tion cannot be calculated without taking into account the 
backreaction due to the perturbation evolution [l,[l7j]- In- 
stead, the background and the perturbations have to be 
evolved simultaneously. In particular, one averages first 
f3T( v \ as in Eq. (|23|) and inserts this into the background 
part of the Klein-Gordon equation (fT6| . In every step, 
the perturbations enter the background equations via /?; 
and Mi. 



C. Simulation scheme 

The methods developed in the previous sections allow 
for a considerable simplification of the numerical treat- 
ment. Rather than evolving a large number of N-body 
particles and the fields p> and ^ona grid, one now merely 
has to evolve a drastically reduced set of differential equa- 
tions. This becomes possible as soon as a collection of 
stable neutrino lumps has formed (at about z ~ 1). The 
preceding cosmological evolution has to be carried out 
with the comprehensive simulation method of Ref. Q. 
Its final state at z 1 provides a distribution of lumps 
at positions xi, with neutrino number factions rest 
masses Mi, and effective couplings This is the start- 
ing point for the simplified scheme. 

Section UlI Bl collects a set of coupled differential equa- 
tions describing the cosmological evolution. These are 
the equation of motion (|15|) . the background Klein- 
Gordon equation (Ti"6j) with its right-hand side (|23|) and 
the usual Friedmann equations. They involve the aver- 
aged potentials at the lump positions, i. e. {Sip(xi)} and 
{^(xi)} (and their gradients) as given by Eqs. ([24)) and 
([2"6[) . Finally, the mass change is computed according to 

dMi 

dt Hi^i dt 



All these equations have mutual de- 



pendences and can only be solved simultaneously. Cold 
dark matter, if included, has to be treated with standard 
N-body techniques. The influence of neutrino lumps on 
the matter component was studied in d, [Tl|, [l8[ . 

The aforementioned equations are only complete to- 
gether with functional relations Pi(fi,<p) and Mi(fi,p), 
cf. Fig. [31 known at all times. As a first approach, one 
may assume a time-independent relation. This is rea- 
sonable if the lumps are virialized and hence their in- 
ner structure is approximately frozen. We will explore 
the stability of individual lumps in Sec. El Furthermore, 
the dependence on (p may be neglected if the derivative 
d{3 / d<p or the variation of ip are sufficiently small. 

It is not clear whether z Ri 1 is late enough for the 
virialization process to have sufficiently proceeded. A 
hint that the cosmological configuration of the neutrinos 
is stabilizing, however, is given by the evolution of the 



total neutrino energy E/ u ) = J d 3 x \fg^ p v (x) oc p„a 3 , 
shown in Fig. |H For a > 0.45, one observes a transition 




0.35 0.4 

scale factor a 



0.55 



FIG. 4. Stabilization of the energy in neutrinos. The dashed 
line shows the evolution calculated by the background equa- 
tions, the solid line is taken from a full simulation run. 

to a regime with a small constant slope. This would be 
compatible with a small monotonic change of the large- 
scale cosmon field and an effective lump mass depending 
on this field, corresponding to the expectation of approx- 
imate mass freezing within neutrino lumps [19]. This 
may be taken as a hint that the neutrino lump fluid may 
become a reasonable picture for a > 0.45. 



ENERGY-MOMENTUM TENSOR OF 
LUMPS 



IV. 



In Sec. IIIH we had to assume properties of the energy- 
momentum tensor associated with neutrino lumps. The 
derivations will be provided in this section. An important 
result is the integrated amplitude A^ v of a single lump's 
energy- momentum tensor, see Eq. ([8]). The derivation in 
Sec. lIV Al includes the vanishing of the total internal pres- 
sure in stable lumps. Next, in Sec. lIVBl we will consider 
the term f3T u , cf. Eq. (f20|) . which sources the energy- 
momentum exchange between cosmon and neutrinos. 



A. Single lump 

We now study a single neutrino lump described by its 
energy-momentum tensor T^ u including contributions of 
the bound neutrinos and the local cosmon field. The 
lump occupies a volume V; its energy-momentum tensor 
vanishes outside. On scales much larger than the lump 
size, it is useful to consider the amplitude 



(27) 



We first switch to the rest frame of the lump where we 
will show A^ v = —Mi 5q 5®. Let us therefor consider the 
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different components separateiy. Cleariy, A°q = —Mi 
by definition of the rest mass. A' t0 — P l is the total 
momentum and thus vanishes in the rest frame, whereby 
A 1 q = 0. It remains to show A l j = 0. 

We assume that the lump is approximately static, i. e. 
its energy-momentum content in a physical volume is 
conserved, 



d (a 3 T» v ) a 0, 



(28) 



a local cosmon contribution) and the cosmon field tp out- 
side the lumps, which mediates the interaction. This 
differs from the usual split in the neutrinos and the 
cosmon T{*K introduced in Sec. UTAl The total energy- 
momentum content T^ t can thus be expressed in two 
ways, 



1 tot — 1 ( u ) ' 1 (ip) ~ 1 lumps 



rpflU 



(33) 



neglecting the metric perturbations. Together with the 
energy-momentum conservation equation, 



= W x T x 1 =d Q T°,+d l T l ? 



a n 
•3-T°„ 



(29) 



0. 



the staticity condition implies d{F l j 

It is convenient to define the three-vector v = 
(T 1 j , T 2 j , T 3 j ) for a given column j . We have just shown 
diw = 0. From now on, we choose i = j = 1 for sim- 
plicity. The amplitude A 1 \ can then be written as 



A\ 



dx I dy dz v\ = a? 



dx I dS-v, (30) 

.I' 



where S x is the slice of V normal to the x direction. Out- 
side the lump, we extend the area S x to a closed surface. 
We can equally integrate over this closed surface since 
there is no contribution outside the lump. We conclude 
that the integral vanishes since divu = inside the en- 
closed volume. This implies A l \ = 0. The derivation can 
equally be done for arbitrary i and j, whereby A l j = 0. 

In the presence of an external cosmon perturbation 5tp 
sourced by other lumps, the energy-momentum conser- 
vation used in Eq. ([29|l only applies to the full energy- 
momentum tensor T^ t including the contribution due to 
Stp. Spatial variations of tp on the scale of the lump are, 
however, small, such that diT^ ot j w diT l j = divv. If 
the external field dtp varies only slowly, the staticity con- 
dition, Eq. (|2"8"|) . applies to the total energy- momentum 
tensor as well. 

The straightforward generalization of the rest-frame 
result gives the amplitude A^ „ = Mi u^u^/j as antic- 
ipated in Sec. IIII Al The lump, on scales larger than 
its size, is described by a standard one-particle energy- 
momentum tensor 



dTMiu^u u S^(x-xi) (31) 



1 Mi 



(32) 



where we have used u° = dx°/dr and 7 = ^f- 



goou 



The neutrino contribution is completely contained in 
Slumps- The cosmon field splits into tp = tp + Stp\ oc , 
and the contribution of the local perturbation 5tp\ oc is 
attributed to the energy-momentum tensor of the lumps. 
The part of the cosmon energy-momentum tensor not 
depending on the local fluctuation Sip\ oc is 



^d x pd x <p 



V(<p) 



(34) 



which corresponds to the standard form of a scalar-field 
energy- momentum tensor. 

Only the total energy-momentum tensor is conserved 
and we want to investigate the energy-momentum flow 
between the components and T' 1 ^ ps . This will yield 
an effective coupling /3; between the lumps and (p. The 
four-divergence of T/^ is 



x T^ = (V x Vxtp-V'(tp)) 



(35) 



In order to evaluate the right-hand side, we employ the 
equation of motion of the full cosmon field tp inferred 
from Eq. (JJ): 



^ x tp-V'{tp)=PT [v) . 



(36) 



Smoothing this relation at the scale A (cf. Sec. IIII Aj) at 
linear order in 5tp\ oc and inserting into Eq. (|35[) yields 



v, T mA 



(37) 



The right-hand side can be expressed in terms of 
lump properties by making use of the conservation equa- 
tion for the total energy-momentum tensor, VaT^ a — 

-V A T^ A ips . The part V\T& mpe can be analyzed in the 
effective description where lumps are treated as point 
particles. The equation of motion (|15[) implies 



VA^ u A mps «- ]T md»<p. 

lumps I 

Comparison with Eq. (|37p yields 



(38) 



B. Smoothed conservation equation 

In the effective description, the two dynamic compo- 
nents are the collection of lumps (with the neutrino and 



lumps I 

which is the relation used in Sec. IIII Bl 



(39) 
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V. ASPECTS OF STABILITY 

The effective description of the cosmological dynam- 
ics outlined in Sec. IIIII relies on the assumption of stable 
lumps. At the current stage of the comprehensive sim- 
ulation method [8\, however, it is not possible to track 
the evolution of lumps after z 1. In this section, we 
sketch some analytic arguments why stable lumps are 
expected to form. We start with considerations concern- 
ing the angular momentum, Sec. IV Al and construct an 
explicit example of a static configuration using hydrody- 
namic equations in Sec. IV Bl In the following, we neglect 
the metric perturbations. 



when approaching the center, this even amplifies the an- 
gular momentum barrier. 

Of course, these results for a test particle in a central 
potential need not generalize to a distribution of particles 
forming a lump. There, we define a neutrino angular 
momentum density as in special relativity, 



X X J 0)> 



(43) 



which, without the cosmon-neutrino coupling, would sat- 
isfy a conservation equation V Q (Z^"/a) = due to the 

conservation equation for T?y Here, derivatives are 
taken with respect to comoving coordinates. Defining 
the total spatial neutrino angular momentum 



A. Angular momentum 



(44) 



The cosmon-mediated fifth force felt by the neutrinos is 
stronger than but in some respects similar to gravity. The 
field equation for Sip, Eq. p7|) . can be compared to the 
usual gravitational Poisson equation. The cosmon per- 
turbation Sip thus plays the role of a potential - similar 
to the gravitational potential - in which a neutrino parti- 
cle moves. In contrast to the gravitational case, however, 
the particle changes its mass m v = m u {ip) while moving 
with velocity w M according to 



-/3m„ 



u x d\p 



(40) 



The loss of mass when moving towards a minimum of the 
potential implies, by momentum conservation, an addi- 
tional acceleration [8j. Hence, it has to be investigated 
whether neutrino lumps are unstable, i. e. continuously 
shrink to smaller sizes until they are stabilized, e.g., by 
the degeneracy pressure [20l ]. 

The cosmon-mediated fifth force, despite the mass vari- 
ation along a particle trajectory, shares an important 
property with gravity: the conservation of angular mo- 
mentum. For example, a single particle moving in a 
spherically symmetric and static cosmon potential ip(r) 
(in physical coordinates) has the conserved angular mo- 
mentum 



L 



7mj, r 



(41) 



in polar coordinates (r, 9) and with the Lorentz factor 7. 
The equation of motion, written for the radial momentum 
p r = jm^r, then contains an angular momentum barrier, 
which prevents the particle from falling into the center. 
It reads 



Pr 



L 2 /3m v dp 
dr 
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(42) 



This is analogous to Newtonian gravity with an angular 
momentum barrier cx L 2 /r 3 and an inward potential gra- 
dient. The only difference is the variation of m„ (and 7) 
along the particle's trajectory. Since the mass decreases 



the conservation equation for IV? in the uncoupled case 
translates to the conservation law 

d t (a a L«) = 0. (45) 
With the coupling, Eq. ©, we instead obtain 
V Q (a-H 1 ^ ) = -a-^T^ (x l d>p ~ xf&ip) , (46) 
and thus 

d t (o^) =d t J d*x ^ a 2 l%] (47) 

d 3 x a 2 (3T {v) (xWip - x 3 d l p) . 

(48) 

For a spherically symmetric lump and thus cosmon po- 
tential ip = ip(t,r), it is straightforward to show 



x l d J ip — x^d l p = 0. 



(49) 



In this case, the quantity a LVs is indeed conserved. 
This is related to the fact that a spherically symmetric 
scalar field does not carry spatial angular momentum, 



#° = -cp (x*gi<p - aJ'dV) = 0. 



(50) 



The conservation of the total angular momentum a 2 L^ t 
then reduces to the conservation of a 2 L^y 

Our considerations hold for an arbitrary isotropic and 
homogeneous background metric. Thus, a does not need 
to be the cosmic scale factor but can also describe some 
local properties of the metric. Fluctuations of the metric 
around the background metric as well as fluctuations of 
the cosmon around an averaged field as tp in Eq. (|10[) can 
be added to the neutrino energy-momentum tensor in 
Eq. flU} . The right-hand side of Eq. (g5J involves then 
instead of ip and fii instead of /?, resulting in a reduction 
of the change of angular momentum. We conclude that 
angular momentum conservation is similar to standard 
gravity. This constitutes a strong hint for a dynamic 
stabilization of the lump. 
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B. Hydrodynamic balance 

We will now study a neutrino lump within a hydrody- 
namic framework and derive a balance equation for a sim- 
ple class of lumps. For this purpose, we will employ mo- 
ments of the neutrino phase-space distribution function 
f{t,x l ,pj) describing the distribution of particles with 
comoving position x l and momentum pj = m v Uj. A dis- 
cussion of stability based on the Tolman-Oppenheimer- 
Volkoff equation can be found in Ref. [2~i| . For simplicity, 
we will restrict ourselves to first-order relativistic correc- 
tions in this section. The equations of motion for a neu- 
trino particle under the influence of the fifth force can 
then be written as 



V . 



1 



P k Pk 
2m?, 



(3m u djip. 



(51) 



The fully relativistic equation in terms of the four- 
velocity u M is presented in 

We will consider the following moments of the phase- 
space distribution function /: 



n = J d 3 pf(t 7 x 7 p), 
nUi = \ d 3 p-^—f(t,x,p), 



nU l U 3 = / d A p 



Pi Pj 



am v am, v 



f(t,x,p) 



(52) 
(53) 



(54) 



The quantities n(t, x), U(t,x), and <Tij(t, x) are inter- 
preted as the number density, the locally averaged pecu- 
liar velocity, and the velocity dispersion tensor, respec- 
tively. Their evolution equations can be derived from the 
principle of particle conservation in phase-space, which is 
expressed by the continuity equation 



/ + 



d(fx*) , d(fpj) 



dx l 



d Pj 



0. 



(55) 



The whole procedure is similar to the standard case of 
gravity (cf. Ref. [221 ]) with the peculiarity of a varying 
mass m v = m u (ip). 

Integrating over the momentum in Eq. (|55|) yields the 
zeroth moment 



h + d ri (nUi) = 0, 



(56) 



with d Ti = a~ 1 d/dx t . A static number density profile, 
ri = 0, is realized if the microscopic motion adds locally 
up to zero, [7 = 0. This is the case for a locally isotropic 
velocity distribution. For this class of lumps, the equa- 
tion for U , which follows by taking the first moment of 
Eq. ([5"5j) and using the equations of motion (|5"Tj) . takes a 
particularly simple form: 



Ui 



1 



drjCTi] + /3d ri ip 1 



3(7 

2n 



1 



aij/3d rj ip, (57) 



with a = er\/3. 



For a static lump, we demand U — in addition to 
h = 0. A glance at Eq. ([5T]) shows that this requires a 
certain balance between the effective pressure cx d rj (Jij, 
generated by the microscopic neutrino motion, and the 
fifth force oc j3d ri ip. Assuming spherical symmetry, the 
balance equation reads 



V = /v(i-f) 

n V 2nJ 



(58) 



with a prime denoting derivatives with respect to the 
radial coordinate. Here, we have used cry = c&y. Solv- 
ing this equation together with the (radial) Klein-Gordon 
equation for the cosmon field yields static lump configu- 
rations. These lump configurations differ from the solu- 
tions discussed in Ref. [2(| since the stabilizing pressure 
is now provided by the neutrino motion rather than by 
the degeneracy pressure. 

At first sight, it is not clear whether this staticity con- 
dition constitutes a stable equilibrium. We perform an 
exemplary numerical check by simulating a single, iso- 
lated lump with the N-body technique [8|]. Rather than 
starting with a static lump configuration by Eq. (|58l) . 
we use a somewhat smaller velocity dispersion a. In 
the subsequent evolution, the lump shrinks and the neu- 
trino pressure increases. Figure [5] shows how the neutrino 
profile becomes more concentrated and indeed stabilizes. 
The simulation evolves the lump in physical time t. For 
convenience, we have translated time intervals At to scale 
factor intervals Aa by the Hubble parameter at z = 1. 
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FIG. 5. The radial profile oc Anr 2 n(r) of neutrinos inside the 
perturbed lump normalized to unity. 

The pressure cancellation between the contributions 
of neutrinos and the cosmon perturbations, cf. Sec. IIIBl 
is established during the stabilization process. This is 
shown in Fig. [6] Similar to Fig. [2] we observe that the 
pressure cancellation is only established at rather large 
distances from the lump. At smaller distances, a residual 
positive pressure remains. 
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FIG. 6. The total pressure integrated to r = 100 /i _1 Mpc 
(red, solid) and 20/i _1 Mpc (dark red, dashed). For com- 
parison, we also plot the neutrino contribution covering all 
neutrinos (black, dotted). 



VI. CONCLUSION 

We have shown that a simplified, effective description 
of the cosmological dynamics in the growing neutrino 
quintessence model is possible. It bases upon describing 
stable cosmon-neutrino lumps as nonrelativistic particles 
with an effective interaction. After the main idea was 
given (Sec. IIIBI) . several aspects needed to be investi- 
gated. 

The first issue concerns the stability of the lumps. 
We have shown in a hydrodynamic analysis of spheri- 
cally symmetric lumps that the neutrino velocity disper- 
sion indeed stabilizes the lumps against the attractive 
cosmon- mediated fifth force (Sec. lVBj) . On more general 
grounds, stability of the lumps is already expected by an- 
gular momentum conservation which holds similarly to 
the gravitational case (Sec. IV A) . Stable lumps may then 
be characterized by the amount of bound neutrinos. In 
numerical simulations of growing neutrino quintessence, 
we have found lumps containing a fraction up to > 10~ 3 
of all neutrinos in the Hubble volume, reaching a mass 
of ~ 10 17 solar masses (Sec. IIII Al) . The total number of 
identified lumps in the Hubble volume is of order 10 . 

Second, it is not clear a priori that the lumps can be 
described as particles. The most important aspect here is 
the vanishing of the total internal pressure. The neutri- 
nos, however, have reached high velocities and an equa- 
tion of state w v fa 0.1. We have shown that - under 
idealized conditions - the neutrino pressure is exactly 
cancelled by a negative pressure contribution from the 



local cosmon perturbations (Sec. IIV[) . A numerical check 
is given in Fig. [2] Under realistic conditions, the pressure 
cancellation may not hold exactly but to a good approx- 
imation. Approximate cancellation of neutrino and cos- 
mon pressure occurs at a characteristic radius n that is 
substantially larger than the radius of the neutrino core 
of the lump. For an effective particle description, ri is 
the size of the lump. A fluid description requires that 
the typical distance between lumps exceeds n. 

Third and finally, a description of the cosmological dy- 
namics requires the equation of motion for the lumps and 
the field equation for the smoothed field (p mediating the 
interaction between the lumps. These equations have 
been derived in Sec. IIIII The decisive quantity character- 
izing the lump interaction is the effective cosmon-lump 
coupling Pi. For small lumps, it approaches the funda- 
mental coupling P quantifying the cosmon-mediated fifth 
force between neutrinos. For big lumps, the effective cou- 
pling Pi is suppressed by a factor of two to three as com- 
pared to p. Since the attractive force is proportional to 
the squared coupling, this corresponds to a suppression 
of the attraction by one order of magnitude. 

The effective description of growing neutrino 
quintessence complements sophisticated numerical 
techniques as it provides physical insight into the 
dynamics. Furthermore, the effective description could 
prove useful in understanding the evolution for redshift 
z < 1, where numerical simulations have not yet been 
successful [1, [ll| . Quantitative results for low redshifts 
are needed to eventually confront growing neutrino 
quintessence with observational constraints. 

The process of lump formation (in the redshift range 
z « 1 to 2) is very complex and still requires a thor- 
ough numerical treatment. It constitutes, however, only 
a transitional period. Thereafter, a physically sound and 
much simpler picture seems to emerge. 

The concepts and methods of averaging developed 
in this paper are quite general for describing lumps in 
the presence of long-range interactions mediated by a 
field. For our purpose, we needed a relativistic treat- 
ment. A nonrelativistic version may be applied to dif- 
ferent clumping processes where separated interacting 
lumps are forming. An example could be a gas or liq- 
uid of macromolecules for which internal structure does 
not play a decisive role. 
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